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ABSTRACT 

We calculate X-ray properties of present-day galaxy clusters from hydrodynamical 
cosmological simulations of the ACDM cosmology and compare these with recent X- 
ray observations. Results from three simulations are presented, each of which uses the 
same initial conditions: a standard adiabatic, Non-radiative model, a Radiative model 
that includes radiative cooling of the gas, and a Preheating model that also includes 
cooling but in addition impulsively heats the gas prior to cluster formation. At the 
end of the simulations, the global cooled baryon fractions in the latter two runs are 
15 per cent and 0.4 per cent respectively which bracket the recent result from the 
i^T-band luminosity function. We construct cluster catalogues which consist of over 500 
clusters and are complete in mass down to 1.18 X lO^^h-^MQ. While clusters in the 
Non-radiative model behave in accord with the self-similar picture, those of the other 
two models reproduce key aspects of the observed X-ray properties; the core entropy, 
temperature-mass and luminosity-temperature relations are all in good agreement with 
recent observations. This agreement stems primarily from an increase in entropy with 
respect to the Non-radiative clusters. Although the physics affecting the intra-cluster 
medium is very different in the two models, the resulting cluster entropy profiles are 
very similar. 
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1 INTRODUCTION 

Because of numerical limitations, most previous simulations 
of clusters of galaxies have been adiabatic, resulting in ap- 
proximate self-similar scaling for cluster properties. How- 
ever, these models are unphysical in that the cooling time of 
the gas in the central regions of the clusters is less than their 
age; also observed clusters do not scale self-similarly. In this 
paper, we generate mock cluster catalogues from simulations 
that include radiative cooling and we show that these can 
reproduce the observed scaling relations. 

The observed X-ray luminosity-temperature {Lx — Tx) 
relation has a slope as large as 3 (e.g. Edge & Stewart 1991), 
steeper than the predicted value of 2 from self-similarity. 
Cooling flows, which are known to exist in a large popula- 
tion of clusters (Allen & Fabian 1998) , are one cause of this 
as the more massive, cooling flow clusters have an excess of 
luminosity. However, when the cooling flow components are 
corrected for, the observed slopes are still significantly larger 
than 2 (Allen & Fabian 1998; Markevitch 1998). The depar- 
ture from self-similarity is even more prominent in low mass 
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clusters or groups (Ponman et al. 1996; Xue & Wu 2000) 
where the slope can be as steep as 8. These departures from 
self-similarity can help us to learn about the history of the 
intra-cluster medium (ICM). 

The Lx ~ Tx relation predicted from self-similarity, 
when bremsstrahlung is assumed to dominated the emis- 
sion, can be written as Lx oc /gas(l + Ziormf'^'^T'^ where 
/gas is the gas mass fraction and zform is the formation 
redshift. Mohr, Mathiesen & Evrard (1999) and Arnaud & 
Evrard (1999) present evidence that the intracluster medium 
is more extended in low-mass clusters so that /gas is a mildly- 
increasing function of temperature, but the latter also stress 
that this is statistically inconclusive due to the fact that 
their sample is not homogeneous and the results are model- 
dependent. 

One explanation for the above observations is that 
galaxy formation is more efficient in poor groups than clus- 
ters. David et al. (1990) show that the ratio between the gas 
and stellar mass increases with gas temperature. They ar- 
gue that as mergers take place to form massive clusters, the 
gas is heated and this suppresses galaxy formation. Bryan 
(2000) compiled observed mass fractions in groups and clus- 
ters from various studies. By assuming variations of galaxy 
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formation efficiencies and applying it to a simple model, he 
successfully predicted the Lx — Tx relation. This motivates 
our Radiative simulation in which gas in removed from the 
ICM by radiative cooling. We are somewhat fortunate in 
that the finite resolution of the simulation limits the amount 
of cooling (or else the vast majority of the gas would have 
cooled to low temperatures) and leads to a similar amount 
of star formation as Bryan suggests is required to reproduce 
the observations. 

Kaiser (1991) and Evrard & Henry (1991) suggest an 
alternative way of expelling gas from clusters via energy in- 
jection at early times, perhaps feedback from galaxy forma- 
tion. The gas does not have to be removed from the cluster 
entirely but only from the central regions that dominate the 
X-ray emission. Observational evidence is provided by Pon- 
man. Cannon & Navarro (1999) who show that the surface 
brightness profiles of cool and hot systems are not scaled 
version of each other and that the core entropies of low tem- 
perature clusters are higher than can be achieved by gravi- 
tational collapse alone. It appears that the core entropy ap- 
proaches a certain value, approximately 100 /i"^^'' keVcm^, 
at low temperature that they designate the 'entropy floor'. 

For this preheating model to work, it has to happen in 
an optimal way, i.e. the right amount of energy has to be 
injected at the right time. There have been various studies 
which suggest that if an energy is injected well before clus- 
ter collapse, 0.3 keV per particle is required (Lloyd-Davies, 
Ponman & Cannon 2000). If energy injection takes place in- 
ternally, i.e. within collapsed halos, more energy is required, 
1-3 keV per particle (Lowenstein 2000; Wu, Fabian & Nulsen 
2000; Bower et al. 2001). The efficiency of energy injection 
depends on the density of the IGM: the higher the density, 
the more energy is required to achieve the same entropy 
level. These results motivate our Preheating simulation in 
which, in addition to radiative cooling, the gas is preheated 
by raising the specific energy by 1.5 keV at a redshift of 4. 

It has been shown that there is a need of some form 
of feedback energy into the ICM in order to get a suitable 
cooled baryon fraction (White & Frenk 1991, Cole 1991, 
Blanchard, Valls-Gabaud & Mamon 1992). Without such 
process, too much gas would be cooled and converted into 
stars at high redshift. As halos at high redshift are small 
and dense, cooling becomes very efficient. A vast fraction of 
gas is cooled into small systems leaving only a small amount 
of hot gas in the halos to form galaxies at later times, the 
so-called 'cooling catastrophe' (White 1992). The picture of 
this process is clearly in contradiction with observations as 
large amounts of hot baryons are observed in X-ray clusters 
at the present day. 

Unfortunately, an accurate determination of the cooled 
gas fraction in clusters of galaxies is hard to obtain. From re- 
cent observations of the K-hanA luminosity function, Balogh 
et al. (2001) show that the global cooled baryon fraction 
is only about 5 per cent. Using a sample of clusters com- 
piled by Roussel et al. (2000) and Carlberg et al. (1996), 
they estimated a cooled fraction of around 20 per cent for 
systems with kT =lkeV, decreasing to around 10 per cent 
at kT =10keV. Note, however, that these results rely on 
extrapolating the observations out to the virial radii (by as 
much as a factor of 5 in radius for low-temperature systems). 
These observations do not preclude excess cool material in 
a form other than stars, such as neutral or molecular gas. 



Indeed, the models of Bryan (2000) and Wu & Xue (2002b) 
both required a cooled gas fraction in clusters of about twice 
the above in order to reproduce the X-ray scaling relations. 

Despite this uncertainty in the cooled gas fraction in 
real clusters, our Radiative and Preheating simulations are 
likely to bracket the true value. We will show below that in 
all other resepects the clusters in our simulations reproduce 
the X-ray properties of clusters extremely well. This agree- 
ment stems from the entropy profiles which are similar in the 
two simulations and much shallower than in a Non-radiative 
simulation. 

In Section 2, we describe the simulations and explain 
how we construct our simulated cluster catalogues. The 
properties of the ICM are presented in Section 3, and the 
temperature-mass and luminosty-temperature scaling rela- 
tions are discussed in Section 4. Finally we summarise our 
conclusions in Section 5. 



2 METHODOLOGY 
2.1 Simulation details 

Simulation data were generated using a parallel implemen- 
tation of the HYDRA code (Couchman, Thomas & Pearce 
1995; Pearce & Couchman 1997), which uses the adaptive 
particle-particle/particle-mesh (AP'^M) algorithm to calcu- 
late gravitational forces (Couchman 1991) and Smoothed 
Particle Hydrodynamics (SPH) to model hydrodynamical 
forces. Our implementation of SPH is similar to that used 
by and discussed in Thacker & Couchman (2001). 

Results are presented assuming a fiat, low-density cos- 
mology, setting the density parameter fio = 0.35, Hubble 
constant h — 0.7l|^, cosmological constant JIa = A/SHg = 
0.65, baryon density fib = 0.019 cold dark matter 
power spectrum shape parameter F — 0.21 and normal- 
ization (78 = 0.9. The initial density field was realised by 
perturbing 4,096,000 (160"^) particles each of dark matter 
and gas from a regular cubic mesh of comoving length, 
100/i~^Mpc. Thus, dark matter and gas particle masses are 
approximately 2.1 x 10^°/i"^Mq and 2.6 x 10^/i"^Mq re- 
spectively. Each simulation was then evolved to z — 0, typ- 
ically taking around 2000 steps, using 64 processors on the 
Cray T3E at the Edinburgh Parallel Computing Centre. The 
gravitational softening length was fixed at e = 50/i~^kpc 
in comoving co-ordinates (equivalent Plummer value) until 
z = 1, after which it was fixed in physical co-ordinates at 
e = 25/i~'^kpc until z = 0. This choice prohibited unwanted 
relaxation effects (see below). 

In this paper we study 3 simulations that differ only in 
the manner in which the gas particles are cooled and heated: 

Non-radiative: This simulation does not contain any radia- 
tive cooling, and heating occurs solely by adiabatic compres- 
sion and shock-heating. As discussed previously (Thomas et 
al. 2001; Muanwong et al. 2001, hereafter T2001 and M2001, 
respectively) this run does not provide an accurate physical 
description of clusters since gas with short cooling times can- 
not cool, but serves as a useful model for comparison and to 
test numerical effects. 

1 We define h = //q/IOO km s-^Mpc-l 
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Radiative: In this run the gas was aUowed to cool ra- 
diatively, using coohng tables from Sutherland & Do- 
pita (1993). We adopted a global gas metallicity of Z — 
0.3{t/to)ZQ, where t/to is cosmic time in units of the present 
value. Material which had cooled (identified as all gas par- 
ticles with overdensities S > 1000 and temperatures T < 
1.2 X lO^K) was identified on each step, and groups of 13 
particles within a softening length were merged into single 
coUisionless particles (hereafter referred to as 'galaxy frag- 
ments'). Fragments could subsequently accrete more cooled 
gas (within a softening length) but not merge with other 
fragments. We do not discuss the properties of these galaxy 
fragments in this paper. 

Preheating: This is identical to the Radiative run, except 
that all gas particles were impulsively heated by adding 
1.5 keV of thermal energy ai z — 4, corresponding to an 
increase in temperature of approximately 1.2 x lO'^K. 

As we shall see, both the Radiative and Preheating runs 
reproduce the observed cluster scaling relations but contain 
vastly different amounts of cooled gas. 



2.2 Cluster identification 

In this paper, we present results only for clusters at 2 = 0. ^ 
Clusters are identified using an identical procedure to that 
used by M2001, as described below. 

First, we create a minimal-spanning tree of all dark mat- 
ter particles whose density exceeds 317 times the mean dark 
matter density in the box. (For a spherical top-hat model, 
this is equal to the mean density within virialised regions, 
although the precise density used does not matter at this 
stage.) We then define a maximum linking length equal to 
0.5 (317)"^/^ times the mean inter-particle separation and 
use this to split the minimal-spanning tree into clumps of 
particles that serve as potential sites of clusters. The centre 
of each clump is defined as the position of the densest dark 
matter particle, then a sphere is grown around all particles 
until radii are found that enclose average overdensities of 111 
(i.e. 317 r2o), 200, 500, 1000 and 2500 relative to the critical 
density, pcr ~ SHq /8nG. Clumps are written into the clus- 
ter catalogues provided that they contain a mass equivalent 
to at least 500 particles each of gas and dark matter within 
these radii, and that their centres are not located within a 
more massive cluster. 

As we are interested in the X-ray properties of clus- 
ters, we did consider using gas particles rather than dark 
matter particles to define the cluster centre. However, we 
decided against this for two reasons: firstly, in runs with ra- 
diative cooling the densest gas particle was often located in 
the outskirts of the cluster as defined by the dark matter; 
secondly, it made it hard to compare clusters from different 
simulations that had similar dark matter distributions but 
very different gas ones. 

The catalogues are complete in mass down to 1.18 x 
lO"/i"^M0. Our final catalo gues for each simulation have 
of the order of 530, 460, 350, 250 and 150 clusters within 
the overdensities mentioned above, respectively. Although 
we do not report on them in this paper, we have also 

^ For clarity, we choose not to discriminate between groups and 
clusters and use only the latter term. 
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Figure 1. Ratios between the specific energy of the gas and the 
dark matter within the virial radii of clusters in the (a) Non- 
radiative and (b) Radiative simulations. The filled traingles cor- 
respond to clusters that show significant velocity substructure, as 
described in the text. 



created cluster catalogues for a large number of higher 
redshifts, and all the catalogues are available on-line at 
http:/ /virgo. Sussex. ac.uk/. 



2.3 Numerical heating 

When undertaking an A'^-body simulation, it is crucial to use 
an appropriate gravitational softening length (see, for exam- 
ple, T2001). Ours is chosen following the Thomas & Couch- 
man (1992) criterion such that the smallest clusters in all 
simulations presented here have 2-body relaxation times of 
at least five times the age of the Universe. Also, for a metal- 
licity Z — 0, the mass of each dark matter particle is approx- 
imately equal to the critical mass, as estimated by Steinmetz 
& White (1997), above which there would be enough spuri- 
ous numerical heating to suppress cooling within the haloes. 
This is the most conservative case, since Z > when the first 
resolved haloes form in the simulations. 

As a check on the degree of numerical heating in our 
simulations, we plot in Fig. ^ the ratio of the total specific 
energy of the gas, tgas, to that of the dark matter, edm, as 
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a function of mass, A/vir, contained within the virial radius, 
Tvir, (defined as the radius which encloses a mean overde nsit y 
of 111 times the critical density, in accord with Section 2.2). 
The specific energy of the gas consists of the kinetic and 
thermal contributions of hot gas particles (throughout this 
paper we use the term 'hot gas' to mean gas particles with 
temperatures in excess of lO^K) with masses mi, speeds Vi 
(relative to the cluster mean) and temperatures Ti: 



2fjim 



(1) 



where /imn is the mean molecular mass which we take to be 
lO-^^g for a cosmic mix of elements. (Note we assume the ra- 
tio of specific heats, 7 = 5/3, corresponding to a monatomic 
ideal gas.) The specific energy of the dark matter particles 
can similarly be written as 



El 2 



(2) 



In the Non-radiative simulation, the ratio of specific en- 
ergies (top panel of Fig. |^) is independent of cluster mass. 
The outhers that have the highest and lowest values of the 
energy ratio in the figures are subclumps which are falling 
into large neighbouring clusters. We illustrate this by using 
solid triangles to denote clusters where the mean relative ve- 
locity of the gas and dark matter is more than 0.1 times the 
velocity dispersion of the cluster (and in fact one of the open 
squares that lie above the mean relation designates a clus- 
ter that only just failed this cut). When radiative cooling is 
included (lower panel), the gas component is found to have 
a higher specific energy than the dark matter and the ra- 
tio increases with decreasing mass; similar results are found 
for the Preheating simulation. This result shows that radia- 
tive cooling leads to mass-deposition and, paradoxically, a 
heating of the residual ICM through adiabatic compression, 
a result first demonstrated in a cosmological simulation by 
Pearce et al. (2000; although see also Knight & Ponman 
1997). In any event, the difference between the two panels 
shows that physical heating processes overwhelm numerical 
ones in our simulations. 

To further test our results, we show in Fig. ^ average en- 
ergy profiles for the 10 smallest clusters in the Non-radiative 
simulation. The specific energies of the gas and dark mat- 
ter follow each other well at all radii and show no evidence 
for an increase in gas specific energy over that of the dark 
matter near the centres of the clusters. 



3 PROPERTIES OF THE INTRACLUSTER 
MEDIUM 

3.1 Baryon fractions 

In this Section, we look at the relative distribution of 
baryons and dark matter in the simulated clusters. Panel (a) 
in Fig. ^ shows the baryon fraction within the virial radius, 
in units of the global baryon fraction, for the Non-radiative 
simulation. Note that, even for the largest clusters in our 
catalogues, the ICM is more extended than the dark matter. 
This is because, for a given specific energy (see Fig. 0) , the 
gas is better able to support itself in the gravitational poten- 
tial than the dark matter. The effect is more pronounced in 
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Figure 2. Averaged profiles of specific energy of the gas, tgas, 
and the dark matter, edmi and the ratio between the two, for the 
10 least massive clusters in the Non-radiative simulation. The 
units are defined in the usual manner of X-ray observations in 
kT units, such that kT = ^fimne.The vertical line indicates the 
softening length. 



smaller clusters although this may be an artefact of poorer 
resolution in these objects. 

In the Radiative run (Fig. ^ panel b) , the total baryon 
fraction in the clusters is almost unchanged. However, the 
mass fraction of hot gas (i.e. the hot intracluster medium 
that gives rise to X-ray emission) is much reduced, espe- 
cially in the low-mass clusters. The Preheating run (Fig. ^ 
panel c) also reduces the hot gas fraction in the centres of 
clusters, but it does so by a different mechanism. Instead of 
turning hot gas into cooled gas (which term we take to in- 
clude stars and galaxies) , it heats it up and expels it from the 
cluster core. In low-mass clusters, the gas is expelled from 
the cluster altogether resulting in a reduced baryon fraction 
within the virial radius, whereas in high-mass clusters it is 
merely redistributed to larger radii but remains within the 
cluster. 

Neither the Radiative nor the Preheating runs give the 
correct fraction of cooled gas. The global cooled fraction 
in the box is 15 per cent for the former and 0.4 per cent 
for the latter. However, we will show later that the X-ray 
properties of the clusters from the two simulations are very 
similar because they both have similar entropy profiles. Only 
when we look at the outer parts of clusters can the two be 
easily distinguished. 



3.2 The entropy of intracluster gas 

For a given gravitational potential, the spatial distribution 
of the ICM in hydrostatic equilibrium is determined by its 
entropy (plus one normalisation condition such as the pres- 
sure at, or the total mass within, the virial radius). Higher 
entropies correspond to more extended gas distributions 
and, in particular, a high core entropy leads to a reduction 
in the gas density near the centre of the cluster and hence a 
much reduced X-ray luminosity. 

In this paper, we follow the usual practice in observa- 
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Figure 3. Mass fractions of liot (solid circles), cold (open circles) 
gas and the total fractions (squares) within the virial radius of 
clusters in the (a) Non-radiative, (b) Radiative and (c) Preheating 
simulations. 



tional X-ray papers and work not with entropy but with the 
closely-related quantity 



sx{r) 



kTx{r) 
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where Ue is the electron density and Tx is the emission- 
weighted X-ray temperature, 



Tx 



T.imiPiAsoit{Ti, Z)Ti 
T,i'mipiKsott{Ti, Z) 



(4) 



Here nii, pi and Ti are the mass, density and temperature of 
the hot gas particles that contribute to the X-ray emission, 
Z is the metallicity, and Aaoft is the cooling function from 
Raymond & Smith (1977) in the soft band, 0.3-1.5 keV. (For 
a fully-ionized plasma with a cosmic distribution of Helium, 
then Ue ~ 0.88p/mH.) 

In the absence of radiative cooling, the entropy of the 
gas can only increase through shock-heating associated with 
mergers and accretion. Within any given cluster, there will 
be a range of entropies, with a positive entropy gradient 
from the centre outwards (this comes about both because 
shock-heating is more effective for material accreted later 
and because any other distribution is convectively unstable). 
For a self-similar cluster population for which we measure 
sx at some fixed fraction of the virial radius (and hence the 
same value of rie for each) , then sx oc Tx ■ In a real cluster 
population this proportionality will not be exact because 
the profiles are not exactly self-similar (the 'concentration' 
varies with mass — see e.g. Navarro, Frenk & White 1997; 
T2001), but it will nevertheless hold to good approximation. 

However, this self-similarity does not seem to extend 
down to low-mass clusters and groups. Ponman et al. (1999) 
and Lloyd-Davies et al. (2000) found that when measuring 
the entropy at a fiducial radius of 0.1 r^ir, which they termed 
the 'core entropy', hot clusters appeared to follow the pre- 
diction from self-similarity whereas cool clusters lay above 
it. (We note that the term core entropy does not imply that 
clusters show a constant entropy in the innermost regions — 
they do not. It is used in a loose term to define the entropy 
outside the central cooling region but still near the centre of 
the cluster.) Although the observations are not very precise, 
a rough interpretation is that the core entropy flattens off 
and approaches a constant value, termed the 'entropy floor', 
in small systems. Ponman et al. (1996) attributed this effect 
to preheating of the inter-galactic medium before cluster 
formation (and thus motivated the Preheating model in this 
paper) . 

In this Section, we calculate the core entropy of our 
clusters in a similar way as in the above observations. We 
average the value of sx in a spherical annulus centered on 
O.lrvir and of width Ar = 0.02 rvir (there are various other 
ways of doing this averaging but they all give similar results) . 
We note that 0.1 rvir is about twice the softening length for 
the smallest clusters in our catalogues. 

The core entropy of the clusters in the three simulations 
are plotted against their X-ray temperatures in Fig. |^. Also 
shown is an approximate fit to the observed relation from 
Ponman et al. (1999, solid curves) and the prediction from 
the non-radiative numerical simulations of Eke, Navarro & 
Frenk (1998, dashed lines). In the Non-radiative run (panel 
a), the core entropy follows the self-similar relation well, al- 
beit with large scatter. On the other hand, the core entropy 
of clusters in the Radiative and Preheating simulations (pan- 
els b and c, respectively) is in good agreement with the ob- 
servations (the solid curves). In the Preheating simulation, 
the excess entropy relative to the self-similar prediction is 
due to the increase in energy of the gas at z = 4 (given that 
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Figure 4. Core entropy at O.lrvir versus X-ray temperature for 
clusters in the (a) Non-radiative, (b) Radiative, and (c) Preheat- 
ing simulations. The straight dashed line is the self-similar rela- 
tionship from the simulations of Eke et al. (1998). The solid curve 
is an approximate fit to the Ponman et al. (1999) data-points from 
their Fig. 2. 



radiative losses in this simulation are neglegible, the excess 
entropy is conserved). The agreement of the Radiative sim- 
ulation with observations, however, is somewhat fortuitous: 
in the absense of some form of feedback machanism, the 
amount of gas that cools to low temperatures is determined 
only by the resolution of the simulation and the imposed 
metallicity of the gas. 



3.3 Entropy and temperature profiles 

Fig. 1^ shows entropy profiles for high- and low-mass clusters 
from each simulation. The entropy is defined as in equa- 
tion except that we use the mass-weighted rather than 
the emission-weighted temperature. The high-mass cluster is 
the third-largest in the box and was chosen because the two 
most- massive clusters show significant amounts of substruc- 
ture within their virial radii; even for this cluster, however, 
there is an obvious subclump at about 0.7rvir that signifi- 
cantly distorts the entropy profile in the Non-radiative simu- 
lation. The low-mass profiles were constructed by averaging 
the profiles of the 10 lowest-mass clusters in each simulation 
that did not show significant substructure — only about half 
the clusters were deemed acceptable. Thus neither panel is 
representative of a typical cluster in our simulations as most 
have much more substructure than is visible here. 

First note that the entropy profiles of the three runs 
roughly match for the massive cluster at radii greater than 
0.3rvir. The differences are due to different amounts of cool 
gas in subclumps that permeate the outer regions of all clus- 
ters. Within this radius, the Radiative and Preheating runs 
have significantly higher entropy than the Non-radiative one. 
For the Preheating run, this excess entropy arises from the 
energy injection &t z — A and is subsequently almost un- 
changed, except for a slight decrease within 0.04 rvir, where 
radiative cooling has had a small effect. The entropy increase 
in the Radiative simulation, although it follows that in the 
Preheating run closely, comes about through an entirely dif- 
ferent mechanism, namely inflow of high-entropy gas to re- 
place low-entropy gas lost via cooling. That the two agree is 
coincidental except in as much as we have chosen the simu- 
lation parameters such that the two runs have similar X-ray 
properties. 

For the low-mass clusters the behaviour is similar, ex- 
cept that the excess entropy extends out to a larger fraction 
of the virial radius. The Radiative and Non-radiative runs 
have converged by r = 2 rvir since cooling is efficient only in 
the centre of the cluster. In the Preheating run, however, the 
excess entropy is still apparent beyond 2 rvir as the entropy 
was increased for all of the gas. 

For all clusters in the Preheating simulation, the ex- 
cess entropy is much less than expected for gas at the 
mean density of the box at the time of preheating (about 
770/i-^/^keV cm ). This is because the intracluster gas at 
z = was already in overdense structures by 2; = 4. For ex- 
ample, the average density of the gas in the 10 least massive 
clusters is about 20 times the mean gas density at 2 = 4. 
Thus, the additional energy injected into the gas only in- 
creased the entropy by approximately 100 h~^^^ keVcm'^, as 
is evident from Fig. ^. 

Fig. ^ shows mass- weighted temperature profiles for the 
same clusters as in Fig. ^ Emission-weighted temperatures 
are similar to these in the inner parts of clusters but decline 
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Figure 5. The entropy profiles in each of the three simulations of 
(a) the third-most-massive and (b) the 10 least-massive clusters 
that did not show significant substructure. The long dash-dotted 
line indicates the gravitational softening length. Bins wrere chosen 
so as to average over at least 32 particles near the cluster core, 
with more further out. 
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Figure 6. The mass-weighted temperature profiles of (a) the 
third most-massive and (b) the 10 least-massive clusters that did 
not show significant substructure. The long dash-dotted line in- 
dicates the gravitational softening length. Bins were chosen so as 
to average over at least 32 particles near the cluster core, with 
more further out. 



more rapidly in their outer regions dues to the presence of 
dense, cool subclumps. Note that the excess temperature 
in the intracluster medium in the Preheating simulation is 
much less than 1 keV/k, except near the cores of the largest 
clusters. Thus, the gas must have cooled adiabatically, flow- 
ing out of clusters and reducing its density, since the time 
it was heated. By contrast, the gas in the Radiative run has 
flowed inwards, raising its temperature by adiabatic com- 
pression. 

In their inner regions, the clusters are approximately 
isothermal. However, the gas temperature declines rapidly 
at large radii, typically dropping by a factor of 2-3 from 
its peak out to the virial radius, in agreement with results 
from previous simulations (e.g. Frenk et al. 1999). Recent 
observational evidence finds no evidence of departures from 
isothermality of X-ray clusters (Irwin, Bregman & Evrard 
1999; White 2000; see, however, Markevitch et al. 1998) but 
only probe out to about 0.4rvir, just where the temperature 
profiles steepen. If these observations can be extended out 
to larger radii then this would be a strong test of the models. 



3.4 Surface Brightness profiles 

X-ray surface brightness profiles for the third-most-massive 
and an averaged low-mass cluster are shown in Fig. ^. Note, 
once again, the significant substructure that is present even 
in these profiles that have been specially selected to be as 
smooth as possible. This substructure is particularly preva- 
lent in the X-ray surface brightness because of its weighting 
with the square of the gas density. 

The total mass in the intracluster medium within the 
virial radius of the massive cluster in the Preheating and 
Non-radiative simulations is roughly equal, yet the Non- 
radiative simulation has a higher surface brightness at all 
radii — this is because of emission from high density gas in 
subclumps that is heated and expelled in the Preheating run. 

The isothermal-/? model (Cavaliere & Fusco-Femiano 
1976) is commonly used to estimate cluster masses. The 
model assumes a halo to be a spherical, isothermal sphere. 
The surface brightness profile of an isothermal gas in hydro- 
static equilibrium within such a halo is 

Ex(i?) = Sx(0)(l + (i?/i?c)')'^'""'"* (5) 
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Figure 7. The averaged X-ray surface brightness profiles in the 
0.3—1.5 keV band in each of the three simulations of (a) the third- 
most-massive and (b) the 10 least-massive clusters that did not 
show significant substructure. The long dash-dotted line indicates 
the gravitational softening length. Bins were chosen so as to av- 
erage over at least 64 particles near the cluster core, with more 
further out. 



where Rc is the X-ray core radius. 

It is not at all obvious that equation ^ provides a 
good description of the surface brightness profiles in real 
clusters. These often fail to flatten in their inner regions as 
much as is predicted by equation (^), so the central regions 
are then omitted from the fltting and the excess emission 
over the isothermal-/? model is attributed to a cooling flow. 
Furthermore, the outer radius to which one can measure the 
surface brightness accurately is only a small fraction of the 
virial radius and depends upon the size and distance of the 
cluster. 

We now estimate the val ues of /3flt for our simulated 
clusters. Then later, in Section 4.1 , we will test the reliability 



of cluster masses deduced using the isothermal-/? model. An 
accurate comparison with X-ray observations would require 
us to produce mock X-ray data from our simulations, with 
appropriate particle backgrounds, etc., and then to run these 
through the X-ray analysis software. We do not yet have 
the tools to do this so instead adopt a simpler procedure 



that is robust even in the presence of substructure. As seen 
in Fig. 1^, the surface brightness profiles of the simulated 
clusters fluctuate, even well within the virial radius. In order 
to smooth the profiles, we integrate them inwards to obtain a 
cumulative luminosity profile. Using the isothermal-/? model, 
this can be expressed as 

Lx{> R) = Lx(0){l + {R/RcfY^^'^"''', (6) 

where R denotes a 2-dimensional, projected radius. We need 
a statistic that is insensitive to emission from the outer parts 
of clusters (where most of the substructure originates) and 
from a possible central cooling flow. One such is 

Lx{> -Rmid) — Lx{> Rout) 



C 



(7) 



Lx{> Rin)~ Lx{> Rout) ' 

which depends only on the flux between Rin and i?out, where 

-Rin < Rntid < R out . 

In our fits, we fix -Rout = 0.4r2oo which is approxi- 
mately equal to the outer extent of observed X-ray surface 
brightness profiles. We also fix J?in = R^ = Rks where _Rks 
is based on the model of Komatsu & Seljak (2001). They 
determined a relationship between the core radius of the X- 
ray emission and the concentration parameter of the dark 
matter density profile (Navarro et al. 1997), 



Rks - 
where 

c = 6 



0.4r5oo/c, 



lO"/i"^M0 



-1/5 



(8) 



(9) 



(We note that that Komatsu & Seljak inadvertently put the 
factor of 0.4 above into the denominator rather than the 
numerator in Equation ^; the expression as we give it is 
compatible with their Fig. 13.) 

Finally, we take -Rmid = 0.15 r2oo and then solve Equa- 
tions ^ and 1^ for with the results shown in Figure ^. 

The first thing to note is the scatter in the calculated 
values of /3fit. This is because individual surface brightness 
profiles are very noisy due to the presence of substructure. 
Secondly, there is a trend for the mean value of /3fit to de- 
crease with decreasing temperature in the Radiative and 
Preheating simulations with the value of approximately 2/3 
in large clusters. The trend is much less apparent in the 
Non-radiative model where most clusters throughout the 
temperature range have consistently high values of /3fit. Ko- 
matsu & Seljak (2001) propose that such a trend may arise 
from observational bias due to the fact that one observes 
out to a larger fraction of the virial radius in higher-mass 
clusters; that is not the explanation here as we have fixed 
Rout = 0.4-R200 (we tried taking a variable outer radius as 
described by Komatsu & Seljak but found that it made lit- 
tle difference to our results). Rather, the difference in /3flt 
values between the Non-radiative and the other two runs is 
a real one, reflecting the raised entropy and hence the re- 
duced density in the inner regions of low-mass clusters in 
the Radiative and Preheating runs. 

If the data were of sufficiently high quality, one could 
evaluate C at two different choices of -Rmid and then solve 
for both Rc and Pat - In practice, there is a strong correlation 
between these two parameters and the solutions sometimes 
return unacceptable large values of Rc (greater than -Rout) 
and correspondingly large values of /3fit. It is for this reason 
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Figure 8. The surface brightness fitting parameter, f3(n, as a 
function of temperature in the 0.3— l.SkeV band, for clusters 
drawn from the (a) Non-radiative, (b) Radiative and (c) Pre- 
heating simulations. 



that we fix the core radius before fitting. Fortunately, Pat 
is insensitive to the precise choice of core radius: we have 
checked that doubling the value of the core radius, i.e. set- 
ting Rc = 2 Rks , does not alter the results significantly. We 
have also taken the more conventional approach of fitting 
functions of the form of Equation ^ directly to the surface 



brightness profiles. This again gives similar results but with 
larger scatter. This is because the surface brightness fitting 
is more sensitive to the presence of substructure at the outer 
edge of the fit than is the method that we describe above. 

The measurement of the Pat parameter as presented 
above is far removed from the usual observational practice. 
Nevertheless, we shall see in Section that the total mass 
within r2oo when determined using the parameter is in good 
agreement with observations. 

3.5 Cooling flows 

There is currently considerable confusion over the obser- 
vational evidence for cooling fiows in clusters of galaxies. 
Recent results show both a deficit of soft X-ray emission 
lines (e.g. Peterson et al. 2002) but also evidence for low- 
temperature gas at approximately lO^K (Oegerle et al. 2001) 
and even large amounts of molecular gas (Edge 2001). Our 
simulations do not have the resolution to investigate cool- 
ing flows in detail. Here we simply wish to test whether the 
standard method of calculating mass-deposition rates gives 
the correct answer. 

In this Section, we calculate mass-deposition rates using 
two methods: (i) by directly measuring the rate at which gas 
cools out of the hot intracluster medium to temperatures 
below 10^ K and (ii) by measuring the emission from within 
the cooling radius, defined as the radius within which the 
gas has a mean cooling time less than 6 Gyr. 

Fig. ^ illustrates the two measures of the mass deposi- 
tion rate for clusters in the Radiative and Preheating sim- 
ulations. The solid circles show the actual mass deposition 
rate averaged over the last few output times in our simula- 
tions, corresponding to a time interval of just under 1 Gyr. 
Although the mass is measured in h~^MQ, we have plotted 
it in units of /i~^M0yr~^ to aid comparision with the pre- 
dicted rate. The particle resolution limits the measurement 
of the mass-depostion rate to be a multiple of approximately 
3/i~^M0yr~^. Note that the mass-deposition is extrememly 
stochastic in that the mass-deposition rate in the first half 
of this time interval typically differs from that in the second 
half by a factor of more than two. 

The open squares show the predicted mass-deposition 
rate, defined as 



M, 



pred 



5fcrx,rcool/2^mH ' 



(10) 



where Lx.rcooi and Tx.rcooi are the bolometric luminosity 
and mean temperature within the cooling radius within 
which the gas has a mean cooling time of less than 6 Gyr. 
The temperature and luminosity here are estimated from 
emission in the soft band, as described in equations (^ and 
(p^. Although there is a large scatter, the predicted and ac- 
tual mass-deposition rates roughly agree with no evidence 
of a systematic bias of one above the other. Because of un- 
certainties in the definition of the predicted cooling rate, 
one should not make too much of this agreement. Neverthe- 
less, it does eliminate the possibility, within the context of 
the models that we simulate, that incorrect interpretation 
of the X-ray observations has led to a vast overestimate of 
the actual mass-deposition rate in clusters. 
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Figure 9. The mass-dopostion rates for clusters drawn from 
(a) the Radiative and (b) the Preheating simulations. The open 
squares show the predicted mass-deposition rate based on the 
soft-band X-ray emission from within the cooling radius; the solid 
circles show the actual mass-deposition rate averaged over ap- 
proximately 1 Gyr. 



4 SCALING RELATIONS 

4.1 Temperature— Mass relation 

The most direct way to compare the simulated temperature- 
mass^ relation between simulated and observed clusters is 
to use the thermal (mass-weighted) temperature of the gas 
within a small region that is well-observed in X-rays. We 
have already done this and presented the results in a short 
paper (Thomas et al. 2002) , in which we compared the sim- 
ulations described in this paper to observations of 5 re- 
laxed clusters using the Chandra satellite (Allen, Schmidt 
& Fabian 2001). In particular, we compared the normali- 
sation of the temperature-mass relation for matter within 
'"2500 (where is the radius of the sphere that encloses 
a mean density of A times the critical density). The Non- 
radiative simulation agrees with previous simulations of that 
kind (e.g. Mathiesen & Evrard 2001) in predicting temper- 
atures that are too low for a given mass, whereas both the 



^ We use temperature— mass rather than mass— temperature be- 
cause we are complete in mass rather than temperature. 



Radiative and Preheating simulations reproduce the obser- 
vations. Unfortunately, spatially-resolved temperature data 
are as yet available only for a few bright clusters and so the 
overlap in mass between the observed and simulated clusters 
is small; nevertheless, there is no reason to expect that our 
simulated results should not extend up to higher tempera- 
tures. 

A much larger body of data exists for emission-weighted 
temperatures of clusters, with poor spatial resolution. In 
this case some form of modelling is required in order to de- 
rive the mass. The emission-weighted temperature, being 
dominated by the high surface brightness, central regions 
of the clusters, does not change very much with ta, but 
the mass does. Generally, one wants to choose as small a 
radius as possible in order to minimise the extrapolation 
outside the region that is well-observed in X-rays. On the 
other hand, theoretical predictions are for the mass within 
the virial radius (A = 111 in this cosmology). In this paper, 
we compromise and use A — 200 as this is the overdensity 
used in two observational papers that we wish to compare 
with: Horner, Mushotzky & Scharf (1999) and Xu, Jin & Wu 
(2001). A third, Finoguenov, Reiprich & Bohringer (2001) 
uses A — 500 but is easily extrapolated to A = 200 (us- 
ing M oc A~^^^ for the isothermal-/? model at radii much 
greater than the core radius). 

kTx — M200 relations for the clusters are presented in 
Fig We use the cooling table of Raymond & Smith (1977) 
to calculate the X-ray temperature in the soft band (0.3-1.5 
keV), as described in Section ^. The open squares show the 
temperature calculated using all the particles in the clusters. 
However, the presence of cold, high-density gas in the clus- 
ter cores (and also in infalling subclumps) gives emission- 
weighted temperatures that are well below the virial tem- 
perature of the cluster. Accordingly, we also show as filled 
circles, for the Radiative and Preheating runs, the 'cooling- 
flow corrected temperature' obtained by omitting emission 



from within the cooling radius, as defined in Section 3.5 



The change is most important for high-mass clusters, for 
which it has the effect of both tightening the relation and 
bringing it closer to the predicted slope of two-thirds. The 
cooling-flow correction does not work very well for clusters 
in the Non-radiative simulation which have cool, dense gas 
at all radii: for this run, therefore, we show instead the effect 
of omitting all gas particles with cooling times shorter than 
6 Gyr. 

Most of our clusters are smaller than those for which 
X-ray masses have been determined. Hence, to facilitate 
comparison with observations, we fit a power-law to the 
temperature-mass relation for clusters more massive than 
5xlO"/i-^M0. Table| lists the normalisation. A, and slope, 
a of the relation 



kTx 



A{M2oo/3 X lo"/r'M0)'/"keV, 



(11) 



where we have chosen to normalise at a mass scale of 
3 X lO'^*/i~^M0 — towards the upper end of our simulated 
catalogue but the lower-end of most observed ones. 

We have divided the Table into three parts. In the first, 
we list results from previous simulations: EMN96 — Evrard, 
Metzler & Navarro (1996); BN98— Bryan & Norman (1998); 
T2001— Thomas et al. (2001); MEOl— Mathiesen & Evrard 
(2001). These use various cosmologies, but fortunately the 
results do not seem to be very sensitive to this. Much more 
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Figure 10. The X-ray temperature versus mass relation for (a) 
the Non-radiative, (b) the Radiative and (c) the Preheating simu- 
lations. The open squares show the total soft-band X-ray temper- 
ature whereas the filled circles exclude emission from (a) gas with 
short cooling times, or (b,c) gas within the cooling radius. The dif- 
ferent lines are from observations (Horner et al. 1999) using mass 
estimates from galaxy velocity dispersions (dashed line), X-ray 
temperature profiles (dot-dashed line), the isothermal- /3 model 
(dotted line), and emissivity profiles derived by surface bright- 
ness deprojection (triple-dot-dashed line). 



important is the resolution of the simulation. Thus MEOl 
have higher resolution than previous simulations (unfortu- 
nately they do not state their precise mass-resolution in the 
paper) and find a lower normalisation; our current simula- 
tions have a higher resolution again and lower the normali- 
sation still further. The reason why the increasing resolution 
lowers the emission-weighted temperature of the gas is the 
presence of cold, dense gas in subclumps. When we exclude 
gas with a cooling time of less than 6 Gyr from the calcu- 
lation then we find that the emission-weighted temperature 
rises once more to a similar value to that found in the earlier, 
low-resolution simulations. 

The middle section of the Table shows results for the 
clusters described in the current paper for emission in the 
soft-band, 0.3-1.5 keV, both with and without the inclusion 
of gas within the cooling radius. The cooling flow correc- 
tion has little effect on clusters in the Non-radiative simula- 
tion, mainly because of the presence of cool gas in infalling 
subclumps (the effect of removing all gas with short cool- 
ing times is much larger and was shown in the first part 
of the Table). On the other hand, gas with a short cool- 
ing time in the Radiative and Preheating simulations resides 
primarily in the cores of large clusters and its omission does 
significantly raise the emission-weighted temperature. In a 
previous paper, Pearce et al. (2000), it was shown that ra- 
diative cooling raises the temperature of the of the intraclus- 
ter medium and we confirm that result. Unfortunately, this 
is cancelled by the lower temperature obtained by moving 
to higher resolution, so that the net effect is to give tem- 
perature normalisations that are little changed over earlier, 
non-radiative, low-resolution simulations 

In the lower portion of Table ^, we present some ob- 
servational determinations of cluster temperatures from: 
HMF99— Horner, Mushotzky & Scharf (1999); FRBOl— 
Finoguenov, Reiprich & Bohringer (2001); XJWQl— Xu, Jin 
& Wu (2001). It is noticeable that various methods provide 
very different scaling relations, both in normalisation and 
slope. The two methods that provide the best agreement 
with the simulations are those that combine optical velocity 
dispersions either with ASCA temperatures or with surface- 
brightness deprojection of Einstein data to create emissivity 
profiles (White, Jones & Forman 1997) . Unfortunately, these 
are the least reliable as mass estimates from velocity disper- 
sions are prone to projection effects (e.g. van Haarlem, Frenk 
& White 1997). Also, the deprojection method requires a 
large extrapolation from A ~ 2000 out to A = 200. The 
highest normalisation is provided by using resolved temper- 
ature profiles. In principle this should be the most accu- 
rate method but as yet the temperature profiles are poorly- 
determined and a high degree of modelling is required. In 
two papers by Nevalainen, Markevitch & Forman (1999, 
2000), for example, the enclosed gas mass fraction in the 
clusters A 401 and A 3571 can be seen to be steeply rising 
at the virial radius, contrary to expectation. 

The greatest degree of consensus is given by different 
authors using the isothermal-/? model, which is another way 
of saying that different groups measure the same relationship 
between /3fit and temperature. In this model, the mass profile 



kT r 

M{< r) « l.llxl0"/3flt — - ^— ^ 
keV ri + r^ 



ft-iMpc 



/i-iMQ,(12) 
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Table 1. Power-law fits to the observational and simulated 
temperature-mass relations of X-ray clusters: cluster sample; 
slope of relation, a; value of fcT/keV at 3 X lOl'*/i~lM0, A. 



Sample 






A 


EMN96 


soft band 


1.50 


3.6 


BN98 


bolometric 


1.50 


2.9 


T2001 


bolometric 


1.50 


3.6 


MEOl 


broad band, high res. 


1.51 


2.6 


Non-radiative 


bolometric. 


1.50 


2.1 




bolometric, tcool ^ 6 Gyr 


1.51 


3.3 


Non-radiative 


soft-band 


1.69 


2.0 




cooling-flow corrected 


1.69 


2.2 


Radiative 


soft-band 


1.96 


2.8 




cooling-flow corrected 


1.64 


3.4 


Preheating 


soft-band 


1 fi'7 
i-.O 1 


Q O 
O.Z 




cooling-flow corrected 


1.61 


3.5 


HMS99 


velocity dispersions 


1.53 


3.8 




temperature profiles 


1.48 


5.1 




emissivity profiles 


2.06 


3.7 




isothermal-/? model 


1.78 


4.4 


FRBOl 


isothermal-/? model 


1.67 


4.1 




polytropic-/? model, low-T 


1.87 


4.3 




polytropic-/? model, high-T 


1.48 


4.4 


XJWOl 


NFW-model 


1.81 


4.9 




isothermal-/? model 


1.60 


4.3 



which gives a mass within r2oo of 
^13 f Pfit kT 



M; 



200 



! 7.69 X 10' 



/ /3flt fcT V 
V keV ) 



'"200 



h-^MQ. (13) 



In most cases r2oo 3> rc so that the correction term in the 
above equation for the finite core radius is approximately 
unity. 

The most extensive analysis of this kind was performed 
by Finoguenov, Reiprich & Bohringer (2001) for the HI- 
FLUGCS (Highest X-ray Flux Galaxy Cluster Sample) from 
the ROSAT All-Sky Survey. They used /J^t values taken from 
fits to the ROSAT PSPC data and temperatures mostly 
from ASCA. Where information on temperature gradients 
was available they generalised the /3-model to allow a poly- 
tropic Equation of State and this gave very similar results. 
A similar result is found by Xu, Jin & Wu (2001) using a 
smaller sample and by Horner, Mushotzky & Scharf (1999) 
using ASCA data on 38 clusters from Fukazawa (1997). 
The slopes of the observed temperature-mass relations are 
slightly steeper than 1.5, refiecting the fact that the mea- 
sured values of /3flt are a slowly increasing functio n o f mass. 



as we found for the simulated clusters in Section 3.4, 

To test whether the temperature-mass relation using 
the isothermal-/? model is consistent with the observations, 
we estimate the masses of the simulated clusters using equa- 
tion ^ and the values of /3fit from Section 3.4 . The resulting 
scaling relation, shown in Figure |ll|, is fully consistent with 
the observations. We do not want to over-interpr et th is re- 
sult as the method of determining Pat in Section 3.4 is far 



removed from the analysis that is carried out on real X- 
ray data. Ideally, one would create mock observations from 
the simulations and analyse them in the same way but that 
is a complex procedure that is beyond the scope of the 
present paper. Nevertheless, we tentatively conclude that 
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Figure 11. The X-ray temperature versus mass relation for (a) 
the Non-radiative, (b) the Radiative and (c) the Preheating simu- 
lations, where the mass is estimated from the isothermal-/? model. 
The dotted, dashed and dot-dashed lines are observed masses us- 
ing isothermal-/? model from HMS99, FRBOl and XJWOl, respec- 
tively. 
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the isothermal-/? model underestimates cluster masses and 
that there is no disagreement between the masses of simu- 
lated and real clusters. 



4.2 Luminosity— Temperature relation 

We define the bolometric X-ray luminosity, estimated from 
emission in the soft band as 



Lx 



Aboi(rx) 
Asoft(rx) 



mipiAsoit{Ti, Z) 



(14) 



where /xmn = 10~^'*g is the mean molecular mass, Tx is 
the soft-band X-ray temperature as defined in equation 
and Aboi and Asoft are the bolometric and soft-band cooling 
functions, respectively. 

Lx is plotted against Tx both with and without emis- 
sion from within the cooling radius in Figure We have 
restricted the temperature ranges to those for which the cat- 



udged by looking at the 
W\). Also shown on the 



alogues are reasonably complete (as 
upper locus of the points in Figure 
Figure are observed relations from Xue & Wu (2000) who 
provided a compilation of observed X-ray temperatures and 
luminosities from the literature. Their sample was divided 
into 3 subsamples: groups (below 1 keV), clusters (above 1 
keV) and the mixture of the two. Their best fit of each cat- 
egory is shown as the dashed, dot-dashed and dotted lines 
in the Figure. 

The first thing to note is that the luminosities of the 
clusters in the Non-radiative simulation are much greater 
than for observed clusters. These are much reduced, how- 
ever, by the removal of gas with short cooling times. At tem- 
peratures above 1 keV, the uncorrected relation follows the 
self-similar relation Lx oc Tx expected for bremsstrahlung 
radiation. At lower temperatures, the luminosity might be 
expected to exceed the self-similar prediction because of the 
added flux from line emission, but in fact it is reduced. This 
is because the cores of the smaller clusters are relatively less 
well resolved. We do not regard this lack of resolution as im- 
portant because the gas in the core has a short cooling time 
and contributes a negligible amount to the total emission in 
the Radiative and Preheating simulations. 

Both the Radiative and Preheating simulations show 
Lx-Tx relations that lie much closer to the observations. 
The agreement is best for kTx >1 keV, less so at lower tem- 
peratures. This may be because we have not raised the en- 
tropy sufficiently in the cores of these systems. However, 
the observational determination of the X-ray luminosity of 
low-temperature clusters is very hard (see e.g. Helsdon & 
Ponman 2000; Wu & Xue 2002a; Voit et al. 2002) and so 
the lack of agreement is not so serious. For both simula- 
tions, but most especially for the Radiative one, omission of 
gas within the cooling radius vastly reduces the scatter and 
brings the outliers down to the main relation. 



5 CONCLUSIONS 

We have analysed the properties of clusters drawn from three 
A'^-body, hydrodynamical simulations of the ACDM cosmol- 
ogy. Each uses the same initial conditions but varies in its 
treatment of the gas physics: a standard adiabatic (Non- 
radiative) model, a Radiative that includes radiative cooling 
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Figure 12. The bolometric X-ray luminosity vs. temperature as 
estimated from the soft band for (a) the Non-radiative, (b) the 
Radiative and (c) the Preheating simulations. The open squares 
use the total soft-band emission, whereas the filled circles exclude 
emission from (a) gas with short cooling times, or (b,c) gas within 
the cooling radius. The dashed, dot-dashed and dotted lines are 
observed relations from Xue & Wu (2000); the solid line in panel 
(a) shows the self-similar relation Lx oc Tj- . 
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of the gas, and a Preheating model that also includes cool- 
ing but in addition impulsively heats the gas prior to clus- 
ter formation. Each simulation generated over 500 clusters, 
complete in mass down to 1.18 X 10^^/i"^Mo. 

The Non-radiative simulation does not reproduce the 
observations but was used as a test of the simulation pro- 
cedure. The clusters drawn from this simulation show no 
signs of numerical heating and behave self-similarly in their 
properties. 

Both the Radiative and the Preheating simulations re- 
produce three key observational relations: 

• The entropy in the cores of low-mass clusters lies above 
the self-similar relation. The measured value at 0.1 rvir tends 
towards a value of approximately 100 /i~^/'^keVcm'^ at low 
masses, with very large scatter. 

• The luminosity-temperature relations are much reduced 
in normalisation relative to the Non-radiative simulation, 
and lie close to the observed relation above 1 keV once cor- 
rected for cooling- flow emission. At lower temperatures we 
still seem to overpredict the X-ray luminosity, although the 
observational errors are large. 

• We have shown in an earlier paper (Thomas et al. 2002) 
that the temperature-mass relation in the inner parts of 
clusters, within r25oo, agrees with observations. In this pa- 
per we reproduce earlier results that show that simulated 
cluster masses within r2oo are significantly greater than ob- 
served ones for a given cluster temperature. However, we 
show that the use of the isothermal-/? and related models 
can lead to an underestimate of cluster masses and once 
this is taken into account the observations and simulations 
are once again brought into agreement. The implications of 
this for the determination of erg from the observed cluster 
temperature function are the subject of a separate paper, in 
preparation. 

The basic explanation for the agreement in the prop- 
erties of simulated and observed clusters is an increase in 
entropy in the cluster cores, over and above that expected 
in an adiabatic simulation. In the Radiative simulation this 
occurs via the removal of low-entropy gas by radiative cool- 
ing, whereas in the Preheating simulation it comes about 
through the imposed energy increase at high redshift. These 
two mechanisms differ considerably in the amount of cooled 
gas that result: at the end of the simulations the global 
cooled baryon fractions are 15 per cent and 0.4 per cent 
respectively, bracketing the observed value. Thus, while nei- 
ther model is a correct description of clusters, one might 
expect that the true model gives rise to similar entropy pro- 
files. 

We showed that the mass-deposition rate in cooling 
flows, i.e. the amount of gas cooling to low temperatures 
in the cluster cores, is reasonably well-approximated by the 
usual method of dividing the luminosity by the enthalpy 
within the cooling radius within which the mean cooling 
time is equal to 6 Gyr. However, the actual mass-deposition 
rate is highly stochastic and may be driven by the infall of 
high-density subclumps. Higher resolution simulations are 
required to investigate this further. 

In the simulations that we have described in this paper, 
the cooling is limited by the numerical resolution. Future 
simualtions will move to higher resolution and will have to 
include the feedback of energy from supernovae. This will 



act as a form of preheating, thus removing the ad-hoc na- 
ture of the current model, although it may be that other 
heating mechanisms are also required. We fully expect that 
realistic models will emerge that successfully replicate all 
the observed features of the intracluster medium. 
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